Quantum state storage and processing for polarization qubits 
in an inhomogeneously broadened A-type three-level medium 
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We address the propagation of a single plioton pulse with two polarization components, i.e., a 
polarization qubit, in an inhomogeneously broadened "phaseonium" A-type three-level medium. We 
combine some of the non-trivial propagation effects characteristic for this kind of coherently pre- 
pared systems and the controlled reversible inhomogeneous broadening technique to propose several 
quantum information processing applications, such as a protocol for polarization qubit filtering 
and sieving as well as a tunable polarization beam splitter. Moreover, we show that, by imposing 
a spatial variation of the atomic coherence phase, an efficient quantum memory for the incident 
polarization qubit can be also implemented in A-type three-level systems. 

PACS numbers: 03.67.-a,42.50.Ex,42.50.Gy 



I. INTRODUCTION 

The propagation of electromagnetic pulses in multi- 
level media has been widely investigated in the last 
decades. Three level atomic media interacting with two 
optical fields in a A-type configuration have been among 
the most considered systems, leading to the discovery of 
a large variety of phenomena, such as coherent popula- 
tion trapping, electromagnetically induced transparency 
(BIT), and slow light [ij. More recently, the development 
of quantum technologies for quantum information appli- 
cations has triggered a renewed attention on the subject. 
In particular, the propagation of weak (quantum) pulses 
in a A medium in the presence of strong (classical) driving 
fields has been considered in detail, giving rise to propos- 
als and implementations for quantum state storage and 
processing 0. 

The non-absorbing propagation of a pair of pulses in 
coherently prepared A-type media has been deeply in- 
vestigated in relation to pulse matching the dark 
area theorem Q, simultons [1| and adiabatons [§, [lO} . 
In all these investigations, the initial quantum state of 
the atomic population plays an essential role. The term 
"phaseonium" was introduced by Scully to describe 
a coherent pure-state superposition between the ground 
levels of a A-medium and methods to prepare a phaseo- 
nium state have been put forward for ladder [12| and 
A [l^ type systems. An interesting effect that takes 
place in these systems is the loss- free propagation of light 
pulses in otherwise optically thick media. In particular, 
it has been shown that matched pulses, i.e., pulse pairs 
with identical envelopes, can propagate without distor- 
tion if the atoms are prepared in a suitable phaseonium 
state 0]- 

Particularly relevant for the purposes of the present 
study is the propagation of two weak optical pulses 
through an inhomogeneously broadened phaseonium 
medium. For this system, it was shown in Rcf. [5| 
that, under the two-photon resonance condition, a cer- 



tain superposition of the fields — the antisymmetric nor- 
mal mode — does not couple to the coherent atomic 
state, whereas the orthogonal superposition — symmetric 
normal mode — does interact with the atoms and is com- 
pletely absorbed. Both the antisymmetric and symmetric 
modes are determined by the phaseonium state, and can 
therefore be tuned according to the phase of the atomic 
coherences. 

In this article, we consider the combination of the 
above mentioned propagation effects occurring in a A- 
type three level medium interacting with two weak pulses 
together with the controlled reversible inhomogeneous 
broadening (CRIB) [l4| technique, first introduced by 
Moiseev and Kroll [15| for quantum state storage. In the 
original CRIB proposal [iJl , a Doppler broadened atomic 
vapor was considered for storage and retrieval of a sin- 
gle photon, exploiting the time-reversal symmetry of the 
optical Bloch equations. Later on, this seminal idea was 
extended to solid state systems with a twofold advan- 
tage: long-living metastable states allow for longer stor- 
age times and the inho mog eneous broadening can be ar- 
tificially controlled [16l - [l9f . CRIB techniques have been 
originally developed for systems composed of two-level 
atoms with the aim of storing quantum information en- 
coded in time domain — so called time-bin qubits. How- 
ever, many quantum information processes and sources 
of photon states arc based on the polarization degree of 
freedom, showing a high degree of interferometric stabil- 
ity and experimental compactness. Most proposals ad- 
dressing the storage and retrieval of polarization qubits 
have considered the possibility of spatially splitting the 
original qubit and storing each polarization component 
separately [l^l- Only recently, it has been discussed the 
implementation of polarization quantum memories that 
do not rely on the spatial splitting of each polarization 
component of single photons [2l| - |23l | . 

In this work, we will show that the CRIB technique 
offers also the possibility to store and retrieve quantum 
states of polarization qubits in a coherently prepared A- 
type system with ideally perfect efficiency and fidelity. 



Wc notice here that the apphcation of the CRIB tech- 
nique to this kind of systems is not straightforward since, 
as aheady mentioned, part of both polarization compo- 
nents of the single photon will propagate as matched 
pulses and, therefore, it will not be absorbed. Hence, the 
photon quantum information can not be, in general, per- 
fectly stored in a A-type medium. Taking profit of this 
fact, we will show that tunable polarization filters and 
sieves can be devised, as well as generic tunable polar- 
ization splitters. Moreover, we will show that, by forcing 
the field components to propagate without the possib- 
lity of adjusting their amplitudes as matched solitons, an 
efficient quantum memory can also be implemented in 
three-level A media. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the physical system and write down the optical 
Bloch equations that govern its evolution. We obtain 
the solutions for the incident field components in Sec- 
tion |lll] and show that, by applying the CRIB technique, 
the system can be used to implement a tunable polar- 
ization qubit splitter. Next, in Section ITVl we consider 
the case for which the phase of the coherent preparation 
of the medium is position dependent. For this particu- 
lar case we demonstrate that, for a large enough optical 
depth, both field components can be efficiently absorbed 
and retrieved preserving the amplitude and relative phase 
between them. The validity of the assumptions made 
for the analytical calculations arc verified in Section |V] 
by numerically integrating the full optical Bloch equa- 
tions. Finally, we summarize the results of this work and 
present the main conclusions in Section IvTl 



II. THE MODEL 



We consider the physical system sketched in Fig. [TJ A 
single pulse, with central frequency uiq and two orthogo- 
nal polarization components (for definiteness we assume 
left and right circular polarizations) propagates in the z 
direction. The pulse interacts with a medium in which 
an artificial inhomogeneous broadening, much wider than 
the spectral width of the pulse, has been created. This 
can be achieved, for instance, via Stark shift of the atom 
levels by an externally applied electric field. Left and 
right circular components of the field interact with the 
|1) — |3) and |2) — |3) transitions, respectively, defining a 
A-type interaction. The atoms of the medium have been 
initially prepared in a coherent superposition of the two 
ground levels |1) and |2}. We assume that the two-photon 
resonance condition is fulfilled for all atoms within the 
inhomogcneously broadened profile. 

The evolution of a single atom of the medium can be 
described, in the rotating wave and electric dipolc ap- 
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FIG. 1. (Color online) (a) Physical system under investiga- 
tion: a single pulse with central frequency luq and two circu- 
lar polarization components enters a medium with transverse 
inhomogeneous broadening, (b) Scheme of the A-type inter- 
action between the single pulse polarization components and 
the three-level atoms of the phaseonium medium, fiis and 
Q23 are the Rabi frequencies of the corresponding interaction 
and A is the two-photon detuning. 



proximations, by the following equations: 
iai3{z,t)ni^{z,t) + c.c, 
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where an is the populations of level Oy is the atomic 
coherence between levels \i) and Vlij = {^ijEij^ jh 

is the Rabi frequency with Eij being the slowly varying 
electric field amplitude of the light component coupled 
with transition \i) o dij is the dipole moment of the 
corresponding transition, h is the reduced Planck con- 
stant, and a; 



j — uJi is the transition frequency be- 
tween levels \i) and Note that we consider a closed 
atomic system satisfying on + 0-22 + 033 = 1. Moreover, 
for simplicity, we have not included any incoherent decay 
term into Eqs. ([1]), since wc assume that the lifetimes of 
the excited level and of the ground state coherence are 
much larger than the whole duration of the interaction 
process. 

In three-level media it is common to define a new basis 
of the system in terms of the bright \B) and dark \D) 
states [23j which, under the two-photon resonance 
condition, are given by 



\B) = -i-(f]^3|i)+n23|2)) and 



fix 



{n*2,\i)-ni,\2)), 



(2a) 
(2b) 
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where 51t — \j\^vi^ + |^^23|^ is the total Rabi frequency. 
For continuous wave light fields, it has been shown that 
the dynamics is restricted to the transition \B) |3). 
In this case, once the system enters the dark-state due 
to e.g., spontaneous emission from the excited level, it 
becomes trapped. However, for light pulses the dynamics 
is different. Indeed, Kozlov and Eberly [^01 showed that 
the dark-state completely controls the system evolution 
at large propagation distances when the medium is in a 
phaseonium preparation [Tlj . 

In order to solve analytically the evolution equations 
of the system we follow the treatment given in Ref. 
and [23! , where the transition operators and the field am- 
plitudes associated to each dipole transition are split into 
the forward and backward modes (denoted by the super- 
scripts f and b, respectively): 



(z,t)e-''(""'*-'=°^) 



(3a) 



, ..^^iz,t)e-^(-"'+'«^>, (3b) 

where luq is the central frequency of the pulse and ko 
the corresponding wavenumbcr. From now on, we set 
p = 3, II, i> ~ 1,2 and fi ^ v. We assume that the initial 
population is distributed only between levels |1) and |2), 
with the excited level |3) being unpopulated. Since the 
interaction involves weak pulses, one can neglect [T^ the 
temporal variations in the level populations as well as in 
the coherence (T12. The validity of these approximations 
is checked and confirmed in Sec. |V] via numerical analysis 
of the full set of the Maxwell-Bloch Eqs. ^ and ©. 

Using the forward/backward mode decomposition 
( ([3a|) and (j3bp ') and the weak field approximation (out- 
lined just above) Eqs. ([!]) are simplified as follows: 

d 



b.f 

dt ''P 



iz,t,A)^iAal' {z,t,A) + 



{z,t) 



ni/)^l/p {z, t). 



^tJ,u\YtJ.iyJ^''^p\^,''J, (4) 

where A = luq — uj^i = ujq — W32. Notice that we have in- 
cluded an explicit dependence on the phase = 4>ij.~4'v 
between the two ground levels through the coherence a^^ . 
This allows to take into account the general case in which 
the system is initially prepared in an arbitrary coher- 
ent superposition of the ground levels. Note that, from 
Eqs. (|4]) it is clear that for atoms prepared initially in 



an incoherent mixture of the ground states, |(Tp 



0, 



the equations for each transition of the lambda system 
become decoupled. This decoupling also occurs for a po- 
larization qubit memory using V-t vpc three level atoms 
as it was recently reported in Ref. [23| . 

The propagation of the forward and backward modes 
of the light, in a reference frame moving with the pulses, 
can be described by the following equations: 

d_ 

dz 



T"7pp 



dA X 



where + (— ) refers to backward (forward) modes, 
Gpp (A) is the inhomogeneous frequency distribution of 
atoms of the corresponding transition; r^^^p ~ g'^Ndf^p/Hc 
with g"^ = uJo/2eoV being the coupling constant; Eq is the 
vacuum electric permittivity, N is the number of atoms in 
the quantization volume V, and c is the speed of light in 
vacuum. It is worth noting that, although the treatment 
is performed in the semiclassical formalism, the linearity 
of the equations ensure the validity of this model also at 
the single photon level. Thus, in what follows, both the 
atomic coherences and the fields could be interpreted as 
classical amplitudes as well as quantum operators. 

From (II])-®, we note that the time reversed (t —t) 
equations for the forward propagating modes are equal 
to the backward ones under a sign change in the detun- 
ings and in the field amplitudes. This symmetry in the 
optical-Bloch equations is indeed the basis for the CRIB 
method. Thus, once the forward propagating input light 
pulse has been completely absorbed, in order to retrieve 
the pulse as a time reversed copy of itself, we need to re- 
verse the detuning (A — > —A). This reversing operation 
can be achieved, for instance, by changing the polarity of 
the device that creates the inhomogeneous broadening. 
At the same time, one has to apply a position dependent 
phase matching operation to transform the forward com- 
ponents of the atomic excitations into backward compo- 
nents, so that the retrieved pulse propagates in the back- 
ward direction. This phase change can be performed by 
transfering the atomic coherences back and forth to an 
auxiliary metastable llil [26j state which also allows for 
longer storage times j27| . 



III. TUNABLE POLARIZATION QUBIT 
SPLITTER 

In the following we will investigate how the three-level 
system discussed in the previous section can act as a 
tunable beam splitter for polarization qubits, a device 
of clear interest in quantum information processing. In 
particular, it will be shown that the polarization basis 
that determines the beam splitter action is defined by 
the phaseonium state. 



A. Quantum filter 

Let us first consider the propagation of the incident 
pulse in the forward direction. For a medium of A-type 
atoms prepared initially in a coherent superposition of 



the ground levels |1) and |2), i.e., a phaseonium state 11|, 
Eqs. (HI- ([5]) can be analytically solved following Ref. Il7 |. 
Thus, we insert the solution of Eqs. (|3]) into Eqs. ([5]), and 
Fourier transform the result: 



d 



[Gp,iA)al'^iz,t,A) 



(5) 



-jy^p^M"- i<l^^^^)Hp^ (^)^^"p {z,uj), (6) 
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where we have used the initial conditions (Tpp(z,i = 
— oo) = and we have defined 

/oo />oo 
Gp,, (A) / e^'-^e'^^drdA. (7) 
-oo Jo 

An analytical compact solution of Eqs. ([S]) can be given 
assuming symmetric transitions, i.e., equal electric dipolc 
moments {\di^\ = |c?23| ~ d) and equal inhomogeneous 
distributions (G31 (A) = G32 (A) = G(A)) yielding 
Vis = '?23 = V and i?32(a;) = H3i(w) = H{uj), which 
reads: 



+ f2i;;(0,c.) e'"^- (e^^^")^ - l) , (8) 



where a{uj) = r]H{ui) is the absorption coefficient. These 
equations describe, for each frequency of the pulse, the 
propagation of the two polarization components of the 
incoming field along the medium. Each solution depends 
on both initial polarization components of the weak pulse 
0^3(2 = 0,aj) and ^^23,{z ~ 0,a;)), which implies that, in 
general, the information carried by each component is 
mixed as the pulse propagates through the A-medium, 
due to the two-photon coherence cri2. 

Moreover, it is easy to see from Eq.Q that, in general, 
the fields are not completely absorbed, even for large op- 
tical depths a{ijj)z. In particular, the two polarization 
components of the pulse can only be completely absorbed 
at large optical depths if the condition 



f]i"3(0,C.) 



Cl 
C2 



(9) 



with c^c* 



fj^^ is fulfilled. However, for arbitrary input 
pulses both components will change their amplitudes and 
phases in such a way that the dark-state becomes pop- 
ulated as they propagate . If the optical depth of the 
medium is large enough, eventually the dark-state will be 
fully populated, so the field components will propagate 
as matched solitons, without any further absorption. In- 
deed, it can be seen from Eq. ^ that the field amplitudes 
will evolve until the relation 



(10) 



is satisfied. Note that relations (jOj) and (|T0|) determine 
the symmetric and antisymmetric normal modes, respec- 
tively, introduced in [sj. In fact, the time reversal sym- 
metry of Eqs. (|4|) and ^ is not sufficient alone to guar- 
antee full recovery of the original input state. Clearly, 
the unabsorbed component of the input pulse cannot 
be recovered back. Specifically, the field that exits the 
medium ends up in a definite polarization state, which 
depends only on the phaseonium preparation. Assuming 
a large optical depth aL — 00, the normalized intensity 



of each of the polarization components at the output of 
the medium z = L, from Eq.([8]), reads: 



(11) 



where ri™(L,w) is the total Rabi frequency of the out- 
put field, as defined in ([2]). This result perfectly agrees 
with (2^, where it was shown that, at the output of the 
phaseonium medium, the amplitude of the field coupled 
with an optical transition depends on the ground-state 
amplitude of the opposite transition. Moreover, one re- 
alizes that in this situation the relative phase between 
the output field components 4>^p{L) — (p^p{L)^ using the 



definition fl™ 

p.p 



of 



, is simply the phase 

the (712 coherence. Considering the input pulse to be a 
single photon in the generic polarization state 

\^^n) =aL\L)+aR\R) , (12) 

entering a medium prepared in a coherent superposition 

|^/'at)=Cl|l)+C2|2), (13) 

we will find at the output the state: 

|'/'L) = (C2-C1)*(C2|L)-Cl|i?)). (14) 

with probability 



PL 



|(V-'m| V'Lt)| 



((Til - 0-12) |a_R|^ + (0-22 - 0-21) la^l^ . (15) 



Note that since part of the incident light pulse has been 
absorbed by the medium, the state (|14p is not normal- 
ized. A convenient picture of the process is that the 
preparation of the phaseonium medium fixes the basis for 
which the incident field is filtered: after the dark-state 
is filled, only the so called antisymmetric normal mode 
(fT4|) propagates without absorption. In turn the latter 
depends only on the phaseonium state, which ideally can 
be tuned at will. In other words, the medium acts as a 
tunable quantum filter. 



B. Quantum sieve 

As the pulse components propagate through the 
phaseonium medium, they adjust themselves to fulfill 
condition (jlOp , which allows absorption- free propagation. 
However, before the dark-state is completely populated, 
part of the field (the symmetric normal mode, given by 
(l9|)) is absorbed by the medium. This information stored 
in the optical coherences can be retrieved back by us- 
ing the CRIB technique. In this subsection we study the 
propagation of the retrieved light pulse in the backward 
direction, which is caused by the sign change of the de- 
tunings and the phase matching operation, performed at 
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time t = 0, once the non-absorbed part of the field has 
left the medium [l7l |. In this situation, the equations for 
the optical coherences and the Rabi frequencies of the 
corresponding field components are: 



)n':,Jz,t)+ia^^^tJz,t) 



(16a) 



dA X 



[Gp^(-A)a^,(z,t,-A)], (16b) 



again with p = 3, /i, = 1, 2 and \x ^ v. The initial and 
boundary conditions at the time of switch (t = 0) are: 



f]^,(L,^) = 0. 



h^,W:^{z,s)\ds, (17a) 
(17b) 

I 



The first two initial conditions (|17ap are obtained from 
the solution of the atomic coherences Q at time t = 0, 
whereas the boundary condition (jl7b[) derives from the 
assumption that, at the time when the detuning is re- 
versed, the non-absorbed field has left the medium. The 
above equations can be solved, as for the absorption 
stage, inserting the solution of (|16ap into (|16bp . and 
Fourier transforming the result. By introducing 



+ 00 



Gpp (-A) 



drdA, 



(18a) 



+ 00 



Gpp (-A) 



^drdA, 



(18b) 

the equations for the backward fields associated with the 
two transitions are: 



+ VppJppi^) 



'pv\ 



.,,)nl(z,^) + a^^nl(z,^) 



,_.,)^:Jz,-uj) + a^^n^^Jz,-uj) 



(19) 



As for the case of absorption, an analytic solution for the above equations can be found. In particular, a compact 
expression can be derived assuming symmetric transitions and inhomogcneous distributions, as done in the derivation 
of (El): 



^lp{^,^) = 



■qj{uj) 



il{F{oj)+H{-u)) 



'pu\ 



(20) 



We note that, similar to the solution for the forward prop- 
agating modes ([5]) , the backward recmited components of 
the field arc also a combination of the initial amplitudes. 
From ((20|) , it is easy to obtain that for a large enough op- 
tical depth the normalized intensity at the output surface 
of the medium (z = 0) reads: 



4(0,^) 



f^up(0,^) 



17^(0, a;)2 



(21) 



where 51^(0, w) is the total Rabi frequency for the back- 
ward components at the output. In the last expres- 
sion, we have considered that the spectral bandwidth 
of the pulse is smaller than the inhomogcneous broad- 
ening, F{lo) ~ H{—io) ~ J{lu)/2 ~ 7r/2, and that the 
absorption coefficient reduces to a ~ i]t: /2. In this sit- 
uation, the relative phase between the field components 
corresponds to the phase difference between the ground 
levels, (/)^p(0) — (/),yp(0) = (p^^^. For the single-photon case, 
the backward retrieved components exit the medium in 
a state 



out I 



(ci+C2)(cJ |L)+cJ|i?)), (22) 



which is orthogonal to (|T4|) and corresponds to the so 
called symmetric normal mode, with probability 



Pi 



(23) 



where P^.^^^ is given in (fTS]) . Therefore, the system acts as 
a tunable quantum sieve that can reemit the sieved state 
on demand. Note that since part of the field is absorbed 
by the medium, the sum of the norms of states (|14p and 
([22|l must be one. 

An example of the sieve protocol is shown in Fig. [21 
where the normalized intensity of each polarization com- 
ponent l^lial^ (black lines) and 1^2231^ (gray lines) is plot- 
ted as a function of the optical distance az. The solid 
lines correspond to the propagation of the polarization 
components of the incident field, Eq. ([8]), whereas the 
dashed lines represent the components of the backward 
retrieved single photon, Eq. ([20l) . In this particular case, 
a medium of optical length aL = 10 has been chosen with 
the atomic population prepared initially in a coherent su- 
perposition fulfilling (Til = 0-6, (722 = 0.4 and (/)i2 = 7r/3. 
The two components of the input light pulse have been 
assumed to have the relative phase of ^13(0) — 023 (0) = 
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FIG. 2. Normalized intensities of the single pulse polarization 
components (see text for definition) coupled with the |1) — |3) 
(black lines) and |2) — |3) (gray lines) optical transitions as a 
function of the optical distance az. The solid lines correspond 
to the propagation of the incident field while the dashed lines 
correspond to the backward retrieved components. The initial 
prepared state of the phaseonium medium is crn = 0.6, (T22 = 
0.4, 012 = 7r/3, whereas the field components are initially 
weighed as 43(0, t^) = 0.9 and r^{Q,tc) = 0.1 with relative 
phase between them being 0i3(O) — 023(0) — 0. 

and initial weighs I'T^iQ^tc) = 0.9 and ^3(0,^^) = 0.1, 
where tc denotes the temporal position of the pulse peak 
in the moving frame. In Fig. [21 we observe how the in- 
cident components change their amplitudes along prop- 
agation and that beyond a certain optical length they 
propagate as matched solitons. The total intensity at 
the output is given by Eq. p5|l . whereas Eqs. pT|) deter- 
mine which fraction of the total transmited intensity is 
in each mode. Since during the transient regime part of 
the field has been absorbed, the CRIB technique allows 
to recover (dashed lines) this information stored in the 
optical coherences. Clearly, since part of the field has 
left the medium, the retrieved intensity is smaller than 
the initial one. In this case, the total retrieved intensity 
is given by Eq. ((23)) while the distribution of the atomic 
populations determine which is the fraction of the output 
intensity associated to each polarization component, as 
given by Eqs. (PT|) . 

To summarize, applying the CRIB technique the 
phaseonium medium can act as a tunable polarization 



qubit splitter: part of the field ends up in the antisym- 
metric (unabsorbed) mode, which exits in the forward di- 
rection, while the rest is absorbed. Moreover, the stored 
antisymmetric mode can be recovered on-demand in the 
backward direction with, in the ideal case, arbitrary time 
delay. This system acts as a state filter or qubit prcpara- 
tor, since one could properly adjust the phaseonium state 
to filter out the non desired polarization components of 
a particular photon state. The same idea works for the 
implementation of a quantum sieve, since only a partic- 
ular desired superposition of the two field components is 
stored and the rest exits the medium. 

IV. QUANTUM MEMORY IN A 
LONGITUDINAL PHASEONIUM 

The analysis performed in the previous section is re- 
ferred to a system with all the atoms prepared in the 
same coherent superposition of their ground levels. We 
have seen that the propagation of a pair of pulses through 
this spatially homogeneous phaseonium leads to an au- 
tomatic adjusting of the pulse components to fulfill con- 
dition (jlOp . This prevents complete absorption and, in 
turn, the straightforward implementation of a quantum 
memory. In this Section we will see that the latter can 
nevertheless be overcome by imposing that the prepara- 
tion of the A-system depends on the position. In par- 
ticular, we consider that the phase of the two-photon 
coherence 0i2 varies with the position along the light 
propagation direction: 

c^i2{z)=e^, (24) 

with 6 being the phase imposed in the medium aX z — L. 
This preparation could be implemented, for instance, by 
applying a linear magnetic field gradient coupling the 
two ground levels such that it produces an opposite Zee- 
man shift of their energies that depends on the magnetic 
field strength. After a certain time the magnetic field is 
switched off yielding a spatial linear phase between the 
two ground levels. If so, the equations describing the evo- 
lution of the system are the same as (l4])-([5]), but now with 
the corresponding position dependence of the two-photon 
coherence phase. Following the same procedure as in the 
previous section, an analytical solution for the field com- 
ponents propagating in a position dependent phaseonium 
medium can be obtained: 



=,(-1)" 



i];;;,(o,c.) cosh 



K+z 
2L 



iO , ( K+z 



-n 



vp 



a{uj)L . 
(0,a;)— — smh 



'K+z' 
2L . 



(25) 



where = \J (a.(iS)LY' — 0^, and for simplicity we have cays to zero for large enough optical depths Q!(w)i, and 
taken the particular case a\\ = (T22 = 0.5. Analyzing this only for = the field components will not be absorbed, 
expression carefully, one realizes that, in general, it de- The physical reason for this result is that the field compo- 
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ncnts arc not able to adjust their amplitudes and phases 
to fulfill condition (fTOjl as they propagate, since the phase 
of the optical coherence between the two lower levels is 
different at each position of the medium. Therefore, the 
dark-state is not populated and consequently each com- 
ponent of the field is absorbed. In particular, if the vari- 
ation of the phase along the medium is large compared 
with the optical depth, i.e., ^ a{uj)L, an exponential 



decay of each component is observed: 

n^-^{z,io):^e-^n%iO,Lo). (26) 

As in the standard CRIB approach, the stored infor- 
mation can be retrieved in the backward direction by 
reversing the sign of the detuning and applying a phase 
matching operation. This leads to a backward propagat- 
ing solution of the field components of the form: 

f2^,(z,c.) = c,,hl^p{o, -Lo) + c^.r2'j;(o, -u), (27) 



where 



~7]J{lu) 



zvF{lj) I.T,(F(i^) + ff (-^)) 

e 2 e 2 



QK.r]{F{uj) + H{-uj)) 
[{9^ + L^jfF{uj)H{-uj)) sinh (W) - {-lYiOQ cosh {W)] sinh ^ 



[(-l)''i6lsinh(VK) + Q cosh (H/)] cosh 



zr,H{-uj) 

e 2 Q K- cosh 
QK-i]{F{uj) + H{-Lu)) 



2L 



- (-l)''ie'sinh 



2L 



e 2 e 2 



K^riF{uj) sinh (W) cosh ( 



K- 



[QriH{-uj) cosh (W) + {-lYi0T]{F{uj) + H{-uj)) sinh (W)] sinh ( — 



e 2 QriH{—uj) sinh 



2L 



(28a) 



(28b) 



and K_ = ^{riH{-uj)LY -9^, W = Q{1- z/L)/2 with 
Q = ^ {7]F{u!)L)'^ — 6^. These cumbersome expressions 
can be written in a much simpler form using the assump- 
tion F{uj) ~ H{—u}) ~ J(cj)/2 ~ 7r/2, which leads to 
K- ~ ~ Q. Moreover, one can realize that for the 
case of optically thick medium {aL oo), i.e. when 
the incident field is completely absorbed, solution ([27]) 
simplifies to 

hl^iz,u) = -n'^^{z,Lo). (29) 

Thus, each component of the field is perfectly retrieved 
with a global phase change of tt. 

Fig. [3] shows an example of the absorbed and re- 
trieved normalized intensity of each polarization compo- 
nent Ifiisl^ (black line) and |f223|^ (gray line) as a func- 
tion of the optical distance az. This plot clearly shows 
that the incident polarization components (solid lines), 
given by Eqs. (P5|) . are completely absorbed at larger op- 
tical distance az. The parameter values considered in 
Fig. [3] are an = 022 = 0.5, 9 = Stt, 0i3(O) - (^23(0) = 0, 
/i'i(0,ic) = 0.9 and ^3(0,^^) = 0.1. The dashed lines 
correspond to the components of the backward retrieved 
components, Eq. (|27p . and we observe practically no dif- 
ference with the incident components. At z = 0, the 
relative phase between the two polarization components, 
not shown in the figure, is the same for the stored (for- 



ward input) and retrieved (backward output) photons. 
Thus, the system is able to store and retrieve on-demand 
a single polarization qubit with unit fidelity. 



V. NUMERICAL ANALYSIS 

The results reported in the previous sections are based 
on the assumption that the atomic populations and the 
two-photon coherence a 12 do not evolve in time, since we 
have considered the weak field approximation. In order 
to check the validity of the analytical approach, we have 
performed exact numerical simulations by integrating the 
full optical Bloch equations ([T]) and ([5|) with a finite dif- 
ference method, similar to the one described in [2^ and 
assuming Gaussian temporal profiles for the pulse com- 
ponents. 

An example of the numerical results corresponding to 
the case of the longitudinal phaseonium medium is shown 
in Fig. m In this figure, the normalized peak intensity of 
both field components is shown as a function of the op- 
tical length with the same parameter values as in Fig. [31 
Comparing these two figures we conclude that the an- 
alytical results are in complete agreement with the nu- 
merical simulations. Moreover, we have checked that, as 
expected, the relative phase between the incident field 



8 



ID; 



in I 2 




FIG. 3. Normalized intensities of the single pulse polariza- 
tion components (see text for definition) |f2i3|^ (black line) 
and |S723|^ (gray line) as a function of the optical distance 
az. The solid lines correspond to the incident field while 
the dashed lines correspond to the backward retrieved com- 
ponents. Initially A- type three-level atoms of phaseonium 
are prepared in state an = (T22 = 0.5 with phase gradient 
of 6 = Stt. The field components are initially weighed as 
l[l{0,tc) = 0.9, 41(0, tc) = 0.1 and the relative phase be- 
tween them is 4ii3{0) — 4>23{0) = both for the forward imput 
and the backward output. 
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populations and the two-photon coherence do not exhibit 
relevant dynamics. This confirms the validity of the weak 
field approximations performed in the previous sections. 
VI. CONCLUSIONS 

We have studied the propagation of a single photon 
pulse whose two polarization components are coupled 
with the two transitions of a coherently prepared A-type 
three-level medium presenting artificial inhomogeneous 
broadening. The propagation effects that normally ex- 
hibit this kind of systems have been used in combina- 
tion with the CRIB technique to discuss potential quan- 
tum information applications. On the one hand, we have 
proposed the use of the A-type system as a quantum 
filter. This proposal is based on the fact that part of 
the incident pulse, i.e., the antisymmetric normal mode 
uniquely determined by the preparation of the atoms in 
the phaseonium state, propagates without distortion. On 
the other hand, we have shown that the orthogonal com- 
ponent associated with the symmetric normal mode can 
be efficiently and completely absorbed and retrieved in 
the backward direction using the CRIB method. In this 
case, the system can be used to implement a quantum 
sieve or, considering both orthogonal modes, a tunable 
polarization qubit splitter. 

Furthermore, we have seen that, by adding a position- 
dependent phase coherence in the phaseonium medium, 
the field components can not populate the dark state 
allowing for a complete absorption of both field compo- 
nents. Then, by applying the CRIB technique, both com- 
ponents can be recovered on-demand, thus implementing 
a quantum memory for polarization qubits. 

Finally the validity of the analytical approach, which is 
based on the weak field approximation, has been checked 
by numerically integrating the full optical Bloch equa- 
tions. The numerical results obtained are in a very good 
agreement with the analytical solutions. 



QZ 



FIG. 4. Normalized intensities of the single pulse polariza- 
tion components coupled with the |1) — |3) (black symbols) 
and 1 2) — |3) (grey symbols) optical transitions as a function 
of the optical distance az after numerical integration of the 
full optical Bloch equations 11} and ((5]). The crosses corre- 
spond to the propagation of the incident field while the circles 
correspond to the backward retrieved components. The pa- 
rameters are the same as in Fig[3] 
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